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We study the Thermodynamic Bethe Ansatz equations for a one-parameter quantum field theory 
recently introduced by V.A.Fateev. The presence of chemical potentials produces a kink condensate 
that modifies the excitation spectrum. For some combinations of the chemical potentials an addi- 
tional gapless mode appears. Various energy scales emerge in the problem. An effective field theory, 
describing the low energy excitations, is also introduced. 
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I. INTRODUCTION 



o 
o 
o 

CN ' Recently V.A.Fateev introduced a series of two dimensional integrable deformations of affine Toda theories pos- 

r^ , sessing the property of dualitju'El. These models exist in two incarnations - fermionic and bosonic - and the strong 

d ' coupling limit of one incarnation corresponds to the weak coupling limit of the other one. Apart frppft being amusing 

^"5 , examples of integrable field theories, these models may have some interesting physical applicationsu^Ll. 

t ' In this paper we consider in some details one of the simplest models of that family. The Lagrangian density of the 

Cn fermionic version of the theory is given by: 

"> : c(f^ = E N^^m^mV^. - ^ii'si^^sf] + 1 [{d,<t>ir + {d.cj^^r] - (i) 



M2 



2e'3(<^2-'^i) + e-^l3'l>i _^ e^/302 



o ■ ^ 

^H ■ 

^\ ' where g — /3^/(47r + 0^) and we can identify two types of fermions V'1,2 and two coupled Toda chains described by 

0^ • the bosonic fields 4>i,2- 

In the bosonic representation the phonon modes, 0i,2, disappear and the Lagrangian density is : 

&: ^'"^ = \ t 1 + I^f2f;,,xs - ^f^^^^ ^ ^^^^ ^ 7V2(xiXi)(X2X2)] (2) 



s=l 



where the coupling constants of the two models are connected by the duality relation: 

7 = 47r//? (3) 

As we have mentioned above, the weak coupling limit of the bosonic theory, j << I {1 ~ g << 1), correspond to the 
strong coupling limit of the fermionic one, f3 >> 1. 

The bosonic form of the model is particularly interesting because the Euclidean version of the theory can be 
interpreted as Ginzburg-Landau free energy of two coupled superconductors. We shall discuss this application in 
greater details later in the text as well as in a separate publication. 

Models (0),(p|) have U(l)xU(l) symmetry and the corresponding conserved charges are given by 

J 2 J [1 + (-f/2yxsXs\ 

Thus one can introduce two chemical potentials and modify the Hamiltonian: 

H^H- IhQi - h2Q2 (5) 

The U(l)xU(l) symmetry of the Hamiltonian reflects in the symmetry of the particle multiplets. In absence of 
chemical potentials the spectrum of the model consists of fundamental particles carrying quantum numbers of U(l) 
groups and their neutral bound states. The two-body scattering matrix of fundamental particles is given by a tensor 
product of two sine-Gordon 5-matrices multiplied by a CDD-factor responsible for cancellation of double poles: 



sinnf + isin(7r/A) 
where the S{X, 9) is the sohton scattering matrix of the sine-Gordon model with 



(6) 



A = 2 



471 + 7^ 



(7) 



and is given by 



Slt;{\:0)^e^'^^'n6a-a6, 



^bb 



-Sab + (1 - Sab) 



sinh(A6') 

sinh[A(6l - i7r)] 

isin(7rA) 



(8) 



with 



+ Sa-bSba{l - Sab)- , r, .„ . .-, 

smh[A(0 — iTTJj 

S,i9) ^ r ,.''-^-'^ sinhtMA-^ - l)/2] 
Jq lu cosh(7rw/2) sinh(7ra;/2A) 

To derive thermodynamic equations we start from the discrete Bethe Ansatz equations which emerge when the 
system is put in a box of size L with periodic boundary conditions: 



-iMsinhe,L^(^^^ ...a„) = n S{ej - 0fc)e(ai, -an) 



k^] 



E^mY^ cosh 9 J 



(9) 
(10) 



where ^(oi, ...a„) is the wave function of n fundamental particles in the isotopic space. The matrix on the right-hand 
side of Eq.(g) is related to the trace of the monodromy matrix r: 



Y[Siej-9k)= Trr(A = 0,;0i,...0„) 



(11) 



Up to the scalar factor given by the product of CDD-factors the latter matrix is equal to a tensor product of monodromy 
matrices of the sine- Gordon models: 



f(A, 01, ...ft„) = fW(A, 01, .. AO X t^^HA, 01, ...^„) 



(12) 



The eigenvalues of the sine-Gordon monodromy matrix are known, hence it is straightforward to write down (^ in 
the diagonal form: 



-iMsinhSji _ 



k^J 



nsinh A[(0j — Ua) + i7r/2] -p-i- sinh A[(0j — Vb) + i7r/2] 
^^^ sinhA[(0j -Ua)-i7r/2] -^J- sinhA[(e'j - vt) -m/2] 



(13) 



-p-i- sinhA[(0j — Ua) + i7r/2] -p-i- sinhA[(uh — Ua) + Itt] 
^J^ sinhA[(6'j - Ua) ~ m/2] ^^^ sinhA[(ufc - Ua) - Itt] 



(14) 



) -I- i7r/2] -p-r sinh A[(uf, - Va) + m] 



nsinh A[(0j — uaj -r i/i / ^j _ -p-i- 
.^^ sinhA[(0j-Ua)-i7r/2] ~ l}^ sinhA[(w6 - Va) - 



(15) 



where: 



SoiO) = e-^^(«) if "tr"w"°n! (16) 

^^ J-J- smhfl + isinfvra/A) ^ ^ 

a— 1 ^ ' ^ 

The numbers m± are related to the conserved charges Qi,2- As we shaU demonstrate later 

m±-^[n-(Qi±Q2)] (17) 

The simplicity of this formula cannot obscure its remarkable meaning: propagating particles are superpositions of 
states localized on different edges of the stripe. This result is valid for other Fateev's models where the Toda array 
consists of more than two elastic chains. 

II. THERMODYNAMIC BETHE ANSATZ (TBA) EQUATIONS 



In order to construct the Thermodynamic Bcthc Ansatz (TBA) equations from ( |l3|Jl4|Jl5| ) we consider the case 
I — g = l/v, with integer v > 2, where the complex solutions of these equations (the strings) have the simplest 
classification. For most of the results the fact that v is an integer is not important and one can generalize them for 
arbitrary g > 1/2 replacing j/ by (1 — g)~^. 

To save space and make our notations more transparent we shall sometimes denote the kernels in the integral 
equations directly in terms of their Fourier transforms. For example, the convolution of two functions g and / with 
the function g having Fourier transform (7(w) will be written as: 

Xj 

de'g{9 - e')f{e') = g * f{9) = {g{u;)} * / (18) 

Thermodynamic Bethe Ansatz (TBA) equations have been introduced inEJu. The free energy of the system is expressed 
in terms of the excitation energies En{0): 

F/L = -TyMn /^cosh01n[l + e-^"W/^] (19) 

In the specific case (n^),(E|) we have only one bound state with energy Ei, and the fundamental particle with energy 
Eo. The function Ei satisfies the following equation : 

Tln[l + e^i(^)/^] - Gn *T\n[l + e-^^^^)/^] = Ah cosh0 

+Gio*Tln[l + e-^°('')/^], (20) 

and is directly coupled only to Eq, which is determined by the equation: 

Tln(l + c^o{e)/T^ _ ^^ ^ rij^fi + c-^oC')/^] = Af cosh0 - h+ + h- ^21) 

+Gio * rin[l + e-^iW/^] - Ts * ln[(l + e^-i(«)/^)(l + e'^^^)/^)] 
The two masses are related by: 

Ml = 2M sin[7r/2(2 - g)] (22) 

and the Fourier transform of the kernel is: 

^ ^ tanh[^c^/2(t.+ l)]sinh[7r(^-l)c^/2(^+l)] 

K+ = Ksyrn + 2ai * s = IH — — (23) 

" cosh(7rw/2) 

where Kgym is related to the two-particle scattering phase : 

K,yr.M = d{9) - ±,±log{So{e)), (24) 

and its Fourier trasform is: 



sinh7r/2(l — l/A)wcosh(7raj/A) 
cosh(7rw/2) sinh(7rti;/2A) 






The equations for e„ and e^„ are: 

£_,, = /i_i^/2 - s *T\n[l + e'^-c— 2)/T] (26) 

£_„= 5n,u-ih-v/2 + s*T ln[l + c-'^-i/'^] [1 + e^^-^+i/^] (27) 

+ 6^-2,nS * rin[l + e-*^-/^] + (5„,is * Tln[l + e"^"/"^] 

e^^h+iy/2- s*T ln[l + e'*—^)/^] (28) 

£« = 5„,^_i/i+i^/2 + s * rin[l + e""-i/'^][l + e'"+i/^] (29) 

+ J,_2,„s * Tln[l + c-^'^/^] + Snss * rin[l + e"^"/^] 

The Fourier transforms of the remaining kernels have the form: 

s{lj) = [2cosh[7r/2A(l -.g)tj]]"^ 

4sinh(g7rti;/2A) sinh[7r(l — (7)ti;/A] 
^"^"'^ ^ sinh(7rc^/2A) 

cosh[(A - l)7rw/2A] sinh[7rw/2A] 



X 



GwH - 



cosh(7ra;/2) 

2 sinh(57r(jj/2A) sinh[7r(l — g)uj/2\] 
cosh(7ra'/2) 



where A is defined in (0) and can be also expressed as A = 1 + i^ ^ . 

The subsequent analysis will show that the fields h± are linear combinations of the chemical potentials /ii^2 

h±^hi± h2 (30) 

(recall Eq.(0)). 

III. PROPERTIES OF THE GROUND STATE 

In this section we consider the zero temperature limit of the TBA equations. We assume that the fields h± are 
strong enough to make Eq negative in some interval [—B,B\. Then it is obvious from Eq.(|2(j) that -Ei > 0: 

El ^ Ml cosh e-Goi* E'f^ (31) 

From Eqs.( p7| , p9| )we deduce that all ej-„ with n ^ v are also positive: 

e±n = -ai * E^;'^ - A„,^_2S * e^"j + {n - l)h± (32) 



where 



A_(c.) = 2coth»^^"^[("~"^^"^"'"^)^"]^^"^["^^"("'"^)"1 (33) 

sinh vVl 



_ sinh[7r(iy-n)»] _ ttuj{1 - g) 

"""^ ^" sinhi/f] ' 2(2 -.g) 

Substituting Eq.(p2) into Eq.(pq,E8h we get the equations for €±^: 



-±L/ 



smh[{iy - 2)f^] 
sinh i/fl 



At last, substituting Eq.(|3|) into Eq.(|2|) we get: 






= h± 



sinhi.f^'^^o 



£;, 



(+) 



iC 



r(-) 



1 



sym ■ 



E;^ ' =Mcosh9- -{h+ + h 



( sinhfi 
[ sinh vft 



(-)4-.(-)l 



* er' + e 



(34) 



(35) 



Eqs.( p4|p5| ) determine the structure of the ground state. We cannot solve them in the general case and therefore we 
shall consider only two special cases: a) ft.+ << ft-_; b) /i+ k, h^. 

In the first case e_y > and ei, is negative inside an interval [— Q, Q], with Q >> B and Q ^ cx) as /i+ — )■ 0. Since 

ejy is small it is convenient to invert the kernel in Eq.(|34|) to get: 

sinh j/f2 



* e(+) + e(-) = vh+/2 + 



1 i,j.(-) 

2cosh[(i.-l)f7] ' ° 



2 sinh f7 cosh [(i^- l)f7] 

For 1^1 >> B it is useful to approximate the right hand side of this equation by its asymptotics: 

1 



where 



and 



9 u< ^^ol f *4"^^-«^exp(-C|e|) 

2C0Sh[(z^ — l)il\ I TT 



'^.-Ce' 



de'Eo[e')e 



iy + 1 _ 2~ g 
u — 1 g 



(36) 



(37) 



(38) 



(39) 




FIG. 1. Schematic diagram of the ground state energies in presence of chemical potential for h+ « h- and large B. The 
same situation for small B will be presented in section VII 



Analogously one can rewrite Eq.(p5h as 



E, 



(+) 



K^E^^ ' = Af cosh e-h^/2- 



2cosh[(iy- l)n] 



M) 



(40) 



where 



K{n) 



sinh ri cosh[(3r/ — 1)17] 



2 sinh i^fl cosh[(i^ — 1)17] cosh[{iy + 1)Q] 



(41) 



.(+) : 



Since el/ is very small, to first approximation this reduces to: 



E, 



(+) 



K^E^o ^ ^ Mcoahe-h-/2 



(42) 



i-), 



For large B we can approximate Eqs.(|36|) and ([43) by the Wiener-Hopf (WH) onesQ. The solutions for Eq {6) 



A-), 



-(+) 



Ef{9 + B) and eV'{e) = tV\e + B) are 



(+)/ 



e[-\o^) 



h. 



4+^H = 



Gi+){oj) 



' Gi-\u:)Gi+\Q) 
h+ 1 



1 



1 



iuj + 1 iu} + 
e-^^ 1 



G(-)(0)-icj + G(-)(iC)C- 



iw 



(43) 
(44) 



where a^ is defined in Eg. (pq) , and Gi (Gi ) and G'-+-'(G*^ •*) are analytic in the upper (lower) half plane and satisfy 
the conditions: 



G(+)(c.)G(-)(c.) = X(c.) 



(45) 



sinh v9 



2sinhricosh[(z^ — 1)17] 
C is defined by Eq.(|3^) and B is determined by the following relation: 



G(+)(i) G(-)(0) 



(46) 



To study the low temperature thermodynamic in Section VI it is useful to introduce also the ground state densities 
for Eq and ei,, a and S respectively, defined by the following equations: 



a'^^'> + K *(t'^ ) == — cosh( 
27r 



A-) 



sinh i>n 



!](+) = 



1 



2 sinh n cosh[(iy - 1)17] J [2 cosh[(zy - 1)17] 

Solutions of the WH equations for a'^-^9) = a^-\e + B) and t{e) = I](6l + B) are 



*a(-) = S° 



a(-'(w) 



t{u;) 



47rG(-)(w)G(+)(i)(icj + l) 



1 



G(+)(- 



w 



-CB 



G(-)(0)(-iw + 0) 'G(-)(iC)(C-«w) 



In the limit b) both e±i, are positive and we have for Eq: 



S^+^ + Ksy„^ *E^Q ^ = M cosh e-h 



Again using the WH method we find the following solution: 



J-) ^ 1 

'° Gl^\Lu)Gi+\0)itiu + l){iu; + S) 



(47) 
(48) 

(49) 
(50) 

(51) 

(52) 



with Gl {uj)GI (w) = Ksymi^)- Substituting this solution in Eq.(|3J) we obtain the following estimate for the gap 
in e±i,: 



KU - h{M/h)^^-3y'^ 



(53) 



In the limit u'^ -> 0{g -^ l)M^ - VMh. 

A schematic diagram of these results is presented in Fig.l and Fig. 2 for the case a) and b) respectively. 




FIG. 2. Schematic diagram of the ground state energies in presence of chemical potential for /i+ ~ /i_ and large B. 



IV. THE LIMIT OF FREE BOSONS 



To check the validity of TBA equations (|2C|-|29D and establish the connection between the fields h± and the chemical 
potentials /ii_2 we consider the free boson limit of the TBAs. It is realized in the case g ^ 1 or v^^ -^ 0. In this case 
the kernels in Eqs.( p6[]29| ) become proportional to the delta functions 

s{9) ^ \5{e) 

and K+ — > /, Gn = Goi -^ 0, such that all TBA equations become algebraic. The mass of the bound state in this 
limit is Ml = 2M and the Free Energy is equal to: 



TM 
F/L = / d6'cosh6lln 

27r 



1 + e 



-Bi/T 



l + e 



'Eo/T 



The solution of TBA equations is expressed in terms of ^ = M cosh6/T: 



(54) 



(55) 



1 + e^±"/^ = 



sinh(/i-tn + a±)/2T 



where 



sinha±/2T 
smhh±/2T 



smhh±/2T 



= Vl + c-Eo/T 



Substituting this into Eq.(21) we get 

Losh{h+/2T)VT+7^^^ + [cosh^{h+/2T) + s[iih^{h+/2T)e-^«/^]^/'^\ 



(56) 



(57) 



which has the foUowing solution: 



l + e" 



-Eo/T _ 



sinh ^ 



(cosh,^ - cosh(/i+/2r) cosh(/i_/2r))2 - sinh(/i+/2r) sinh(/i_/2T) 
Combining it with (pq) in the Free Energy, we get 



i^/L = - — / ^6*^'/ cosh 61 In(jc-) 
27r J 



(58) 



(59) 



(60) 



T={l-e 



-2e^-2. 



sinh^ ^ 



(cosh^ - cosh{h+/2T) cosh(/i_/2T))2 - sinh^(/i+/2r) smli{h^/2T) 



{1 



1-e 



-|-(h+ + h_)/2Ti 



[1-e 



-^+(h++h-)/2T 



][l-e- 



-«-(/i+-?i-)/2Ti 



[1-e 



-^+(/i+-/i_)/2T 



1}"' 



(61) 



The free boson limit of TBAs ( |20fj29| ) does reproduce the free energy of two non interacting complex bosonic fields. 
This indicates that the TBA equations are correct. Another important result is the confirmation of Eqs.(p^,pO|). 



V. SEMICLASSICAL ANALYSIS OF THE LIMIT OF SMALL 7 

Let us consider the Euchdean form of the Lagrangian (pi). Rescaling the fields, Xs ~^ (7/2)As, it becomes: 



rib) _ ±_ 

T 



Ef 



^"''■4^ + M' {A, Ai + AjAj + 2(A, Ai)(A2Aj)) 



(62) 



Then in the limit 7 — > the quantum fluctuations are suppressed and one can approximate the ground state energy 
as minimum of the functional 



E{As,h,) 



T 



M- 



{AiAi + A2A2 + 2(AiAi)(A2A2)}- V - 



A. A,, 



(63) 



We will not address this problem, already studied by FateevEl, but we will use this approximation to to get more insight 
to the problem and give an intuitive picture of the low energy physics. Let us perform the following transformation 
in the form (|62^ 



As = sinhp^e*''''' 



(64) 



Under this transformation the measure transform in Z?pD(y9sinh(2(^) and the Lagrangian density in the new fields in 
presence of chemical potentials assume the form: 



Z:^^^ = — { Zl [(9p/0s)^ + tanh^ Ps{df,ipsf - h, tanh^ p^dQ^p,] 

^ s=l,2 

+A/~[sinh pi + sinh p2 + 2 sinh pi sinh P2]} 
For hi = h and /i2 = 0, which corresponds to the situation b in section III, we get: 

£W = ^{{d^pif +t!,n\i^ pi[{d^^if - hd^^i] + {d^p2f + 



(65) 



1' 



2 „ , •„i 2 



(66) 



tanh^ p2{di_,ip2) + M (sinh^^ pi + sinh-^ P2 + 2 sinh-^ pi sinlr P2)} 



The term h{dQ(pi) can be absorbed by the shift: 



and 



tpi ^ (pi - {h/2)t (67) 

which generate also an additional term —h? /A. Then the Lagrangian density becomes: 

/:W = A{(a^^^)2 + (a^^2)' + tanhVi(5p^i)'+tanhV2(a^<P2)' + K//(Pi,P2)} (68) 

where: 

K// (Pi ,P2)^Vi{pi) + M^ sinh^ p2 + 2M'^ sinh^ pi sinh^ pa (69) 

yi(pi) = Af2 sinh^ pi - h^/Atanh^ pi (70) 

For /i/2 > M, VhJpi,p2) has a minimum at pi = po and p2 — where exp2po = \/8{h/M). In the vicinity of the 
minimum we have: 

Veff (pi =po + x,p2^y)^ hM{V2y^ + VSx^) (71) 

which gives the masses going like \/hM, in accordance with (p3|). 

Thus when the chemical potential exceeds the threshold the kinks condense and gapless Goldstone mode appear. 
To see this explicitly we write pi = po + a; and keeping only quadratic terms in x obtain the following expression for 
the Lagrangian: 

jrib) ^ ^{{d^xf + {Mh/2)x^ + tanh^ po{d^^if + (72) 

{^^,P2f + tanh^ P2{d^(P2)^ + sinh^ P2h/M} 

Here we can identify a gapless mode (pi with velocity equal to the bare one, a gapful field x with the mass tti^ ~ hM/2 
and an effective integrable field theory described by the Lagrangian density: 

^2 = i(-)2{Af sinhV2 + (9mP2)' +tanhV2(9M'/'2)'} (73) 

2 7 

In the case hi ~ /12 corresponding to the situation a) of Section III, the effective potential acquires the form: 

K'}^(pi,P2) = 1/i(pi) + yi(p2) + 2M2sinh2pisinh2p2 (74) 

Repeating the same analysis one finds two independent gapless modes with the same velocity. Then the semiclassical 
approximation gives a correct qualitative description of the low energy behavior of the system. As a byproduct we 
obtain another confirmation of (|30|). 

VI. LOW ENERGY PHYSICS: MASSLESS MODES AND ENERGY SCALES 

Let us go back to the zero temperature limit described in Section III. As we have shown, in the case hi « /12 there 
are two soft modes: Eq and e^. At \6\ < B the mode e^, closely follows Eq, a,t\9\ » B its behavior is determined 
by the asymptotics (p^) and this mode becomes independent. The temperature scale below which the two modes 
decouple can be determined by the value of of Ei, at 6* ^ S: 

T.ep-e^H-B) (75) 

For low temperatures, T < Tsep, the two modes give independent contribution to the thermodynamics. To calculate 
the low temperature free energy: 

F/L = - — M f dOcosheinil + e"^"^''^^^] (76) 

27r J 

we use the Eqs. (M,B2) which at finite temperature have the form: 



rin(l + e^"'^) - K * rin(l + e"^"/^) = McoshO - ^ - (- r-^ -r\ * Tln(l + e^-/^) (77) 

2 [_ 2 cosh(z^ — l)s2 J 

| 2sinh»cosh(l-.)» }*^^'-(^ + -"^)-^^-^^ + - '^^)^^-'-- (^«) 

where a^ is defined in Eg. (|37|) . To isolate the contributions that vanish at T=0 we rewrite (u^ like: 

^^+^ + if * E^''> + {I -K) * Tln(l + e-l^ol/^) = r.h.s. of Eq.(77) (79) 

where / is the identity operator. Using Eg. (|47|) we can rewrite the free energy as: 

F/L^-^M /d6'coshe'(4"^+Tln(l + e-l^«l/^)) = /o + /,. (80) 

271 J 



where 



fo = -T f dea{e)ln{l + e-\^''\/^) ^ -TT^T^/3Vc (81) 

/. = -^/'^^{ 2cosh(Ll)» } -^-H^)Tln(l + e^.^/-) « -.^TV^V. (82) 



The velocities of the two modes are determined by 

J 

'de 



K = El^l[2-^^'™-iB (83) 



s=± 



,d^AO). 



For large B all the quantities appearing in these equations can be calculated using the Wiener-Hopf solutions obtained 
in Section III via the following relations: 

dEo'-^^/de ^ dsQ^^^/de = lim ±(iw)2e(±)(w) 



1=B 



?^0 ^ — ^oo 



a^^\B) = d''^\Q) ^ lim ±iujd^+\uj) 

and analogous relations for Cj^ and S. The results are: 

Vc = Vs = l (85) 

The fact that the two velocities are the same seems to be peculiar to the model in consideration and not a general 
property of the family of Fateev's models. In general Vg depends on g and on the number of Toda chains. 

Thus we have shown that the presence of chemical potentials introduces additional energy scales in the problem 
with non trivial crossovers. At temperatures T <$^ h, e_„ and e„ (n ^ v) modes decouple from the other modes. The 
same happens for Ei at energies smaller then their gap and in some physical regions there is also an intermediate 
energy scale, as we will see in the next Section. As one lowers the temperature below T^ep the gapless modes decouple 
from each other and the picture of two non-interacting Luttinger liquids emerges. 

VII. WEAK COUPLING LIMIT: NON PERTURBATIVE RESULTS 

Now we wish to focus our attention on the limit 5 — > 1 which corresponds to the weak-coupling limit of the bosonic 
theory. Naively one can imagine to obtain perturbative results. It turns out however that, in the presence of the 
chemical potential, this is not the case, the reason being that the kinks condensate survives and affects the excitation 
spectrum. A similar effect is present in the Klein-Gordon limit of the sine-Gordon modelQ. 

To be specific let us consider the situation a) of Section III, all the energies e„ {n 7^ v) have gaps of the order of 
h- and can be omitted at temperatures T << h^. To obtain information on the excitation spectrum one has to 
study the fundamental equations ( |3l| ) and (p2|). The standard relativistic dispersion law: Ei{9) = Mi cosh6', can be 
modified by the last term in ( pi| ) which emerge in presence of kinks. 

10 



In the limit g -^ I, the kernels simplify K{uj) = Ko{uj)[l + 0((1 - gf)), and Gio{uj) = Go(w)(l + 0((1 - gf)), 



where 



47r2 cosh(7rw(A + 2)/2A) 
"^'^^ " ^'^cosh(7rw/2) sinh(7rtj/A) 
TT, , cosh(7rw(A + 2)/2A) 

^ — (1 — g)'-^ 

A cosh(7raj/2) sinh(7r(jj/A) 



(86) 



and 



Go = (1 - g) 



■nij sinh(7rcL;/2A) 
A cosh(7ru;/2) 



(87) 



Again we can study TBA equations with the kernel -^^0(0;) and Gq{lo) in the limits of large and small B. We focus 
the attention on this second case where the effects we want to consider can be seen more clearly. To study Eq.(|3l|) 
for small B we need the kernel Kq{9) only for small 0. This is equivalent to approximating Eq.(^^ for large uj, as: 



Ko{u)^^{l-g)\u\ 

Approximating also the rhs of Eq. ( |42| ) for small 6 we can rewrite it as: 

KqIEq = Me'^/2 - H 



(88) 



(89) 



where H = /i_/2 — M and ? denote the convolution of support {—B, B) and Kq is the Fourier transform of (88). One 
can absorb the constant part of the kernel in Eq introducing the new function Eq = ^(1 — g)Eo. Equation ( [42|) then 
becomes: 



KoiEo = M9^/2 - H, 
Kq is the singular integral operator corresponding to the kernel Kq{lu) = \u!\ and defined: 



E1 



1.010 



1.000 



jzEoi9')d0' 




(90) 



(91) 



EO 




-0.01 



0.20 
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FIG. 3. Dispersion relation Ei (upper figure) in the presence of kinks for B << 1 and g — 0.02, from the numerical solution 
of Eq.(pl|). The dotted line represents the solution in the absence of kink condensate. In the lower figure we plot the gap less 
mode _Eo. 



Using this form of the kernel, Eq. ( P0[) can be reduced to the canonical form: 

where wf have introduced: V{9) = dEo{9)/d9 and g{9) = tt{M9'^ /2 — H). This equation has the following general 
solutiom: 

y(0) = i-Vs2-02p f ao'^^JZ} (93) 



which gives: 



EoiO)=-^^^^^^iB--Or/^ (94) 



where B^/A = 'j^^ — - = A. Analogously the ground state density ( ^ ) turns out to be: 

MX 



<'^-2^^Hl^)^^^' (95) 



From here we find the relationship between h- and the charge: 



2^(1-5) 
The Fermi energy, Ep = Eq{9) as function of the charge if therefore: 



Q = ¥^^ M 



--S<'-""^ 



^0 

M^ 



3/2 

(97) 



We can now use the above results to solve Eq(pi|). Rewriting the Fourier transform of the kernel Go as: 



2 d (7r/2(l-.9))sinh( 

^ 'd9 sinh2 6'+(7r/2(l-5))2cosh2 9 



Go- zi^ ■ ,2n , ,_/on — zwrzziTT^ (9^) 



it is possible to extract the asymptotic behavior of E 



J MiCosh6' + A//3cos(2sin-i(e'/B)),6l< B , . 

^~\MiCosh9 + Ae-<^,9^B ^^ 

The numerical solution is shown on Fig. 3. One can clearly see that even in the weak coupling limit the dispersion 
law is drastically modified at small 9. 

This analysis is supposed to clarify the subject of different energy scales present in the model (cfr.Fig.4). At very 
small doping, Q <S^ Qi = M{g — 1)"^/'^, Ep is smaller then M and there are two scales besides h: M and Ep itself. 
The mode Ei is completely decoupled at low temperature. As you increase the doping you reach an intermediate 
region, characterized by M(l — g)^^^^ ^ Q <C Q2 = M(l — g)~^^ where Ep ~ M but still much smaller then Ei. At 
very large doping we find a very interesting region. As we have seen in Sec. Ill at large Q (or equivalently at large B) 
Ep ~ Q, while £'i(0) ~ ^/Q. Then there is a regime in which Ep is bigger then the gap of the ii^i-mode, and there 
will be temperatures, T <C Ep for which the i?i-mode becomes soft. 
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T 



2M_ 



FIG. 4. A Schematic diagram describing different energy scales as a function of the charge. We can recognize three distinct 
regions as described in Sec. VII. 

VIII. DISCUSSION 

We have considered in some detail the thermodynamics of one of the integrable field theories recently introduced 
by Fateev. Most of the results are quite general and remain valid also for other models of the family. These models 
are interesting for various reasons. First of all they have peculiar mathematical properties mostly related to a dual 
representation of the theory. In addition they probably can find application to various physical systems. Considering 
the model as a (l+l)-dimensional Quantum Field Theory one can interpret it as two one dimensional conductors 
coupled via phonon interaction. The problem of quasi one dimensional systems with electron-phonon interaction is 
a very important one and cannot be approached with perturbative methods, it is then very important to have exact 
results for some specific model. This interpretation of the model and the possible application to physical systems have 
been discussed in refB, for physical regimes different from the one considered here (i.e. for g < 1/2). On the other side, 
the Euclidean version of the model, in the bosonic representation, describes an effective Landau-Ginzburg theory of 
two coupled superconductors. In the specific model considered in this work, the coupling between the superconductors 
is quite simple. Nevertheless for other models of the Fateev's family, where coupling between the fermionic modes is 
achieved via a higher number of Toda chains, the resulting coupling between the two superconductors is via elastic 
modes and then can be interpreted as two layered superconductors separated by an insulating stratum. In relation 
to the high-Tc (cuprate) superconductors an interesting problem is to study the dependence of the superconducting 
properties on the number of insulating layers between the superconducting planes. We will address this problem in a 
separate publication. 

Let us give a brief summary of our results. As shown in Sec. VII the results are non-perturbative also in the weak 
coupling regime and could only be obtained through exact or non-perturbative methods. This is related to the fact 
that the chemical potentials generate a kink condensate that survives in the weak coupling limit. As we have shown, 
for some combinations of the chemical potentials two gapless modes emerge and various energy scales appears, which 
characterizes crossovers between different regimes. Probably the most interesting region of the phase diagram is 
obtained at large doping where the £^i-mode become soft at temperatures for which the £^o-one is still frozen. It 
deserves further investigation. Another important result is the one contained in Eq.(17). This equation tells us that 
charged particles in the Fateev's model are linear combinations of particles located on the edges, which is important 
for qualitative understanding of the physics of the model. 
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